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ABSTRACT 

We present constraints on the mean matter density, fi m , dark energy density, £1de, 
and the dark energy equation of state parameter, w, using Chandra measurements of 
the X-ray gas mass fraction (/ gas ) in 42 hot (kT > 5keV), X-ray luminous, dynami- 
cally relaxed galaxy clusters spanning the redshift range 0.05 < z < 1.1. Using only 
the /gas data for the six lowest redshift clusters at z < 0.15, for which dark energy 
has a negligible effect on the measurements, we measure f2 m = 0.28 ± 0.06 (68 per 
cent confidence limits, using standard priors on the Hubble Constant, H , and mean 
baryon density, H,^h 2 ). Analyzing the data for all 42 clusters, employing only weak 
CO ■ priors on Ho and flhh 2 , we obtain a similar result on f2 m and a detection of the ef- 

>" fects of dark energy on the distances to the clusters at ~ 99.99 per cent confidence, 

CO , with f^DE = 0.86 ± 0.21 for a non-flat ACDM model. The detection of dark energy is 

CO ■ comparable in significance to recent type la supernovae (SNIa) studies and represents 

strong, independent evidence for cosmic acceleration. Systematic scatter remains un- 
detected in the / gas data, despite a weighted mean statistical scatter in the distance 
' measurements of only ~ 5 per cent. For a flat cosmology with a constant dark energy 

equation of state, we measure Sl m = 0.28 ± 0.06 and w — — 1.14 ± 0.31. Combining the 
/gas data with independent constraints from cosmic microwave background and SNIa 
studies removes the need for priors on and Hq and leads to tighter constraints: 

f2 m = 0.253 ± 0.021 and w = —0.98 ± 0.07 for the same constant— w model. Our 
most general analysis allows the equation of state to evolve with redshift. Marginal- 
izing over possible transition redshifts 0.05 < z t < 1, the combined / gas +CMB+SNIa 
data set constrains the dark energy equation of state at late and early times to be 
wq = —1.05 ±0.29 and w c t = —0.83 ±0.46, respectively, in agreement with the cosmo- 
logical constant paradigm. Relaxing the assumption of flatness weakens the constraints 
on the equation of state by only a factor ~ 2. Our analysis includes conservative al- 
lowances for systematic uncertainties associated with instrument calibration, cluster 
physics, and data modelling. The measured small systematic scatter, tight constraint 
on Q m and powerful constraints on dark energy from the / gas data bode well for future 
dark energy studies using the next generation of powerful X-ray observatories, such 
as Constellation-X. 

Key words: X-rays: galaxies: clusters - cosmological parameters - distance scale - 
cosmology: observations - dark matter - cosmic microwave background 
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1 INTRODUCTION 

The matter content of the largest clusters of galaxies is ex- 
pected to provide an almost fair sample of the matter con- 



tent of the Universe (e.g. White et al. 1993; Eke et al. 1998, 
Frenk et al. 1999). The ratio of baryonic-to-total mass in 
clusters should, therefore, closely match the ratio of the 
cosmological parameters fib/fl m . The baryonic mass con- 



© 0000 RAS 



2 S. W. Allen et al. 



tent of clusters is dominated by the X-ray emitting gas, the 
mass of which exceeds the mass of optically luminous ma- 
terial by a factor ~ 6, with other sources of baryonic mat- 
ter being negligible (Fukugita, Hogan & Peebles 1998; Lin 
& Mohr 2004). The combination of robust measurements 
of the baryonic mass fraction in clusters from X-ray obser- 
vations together with a determination of Slh from cosmic 
microwave background (CMB) data or big-bang nucleosyn- 
thesis calculations and a constraint on the Hubble constant, 
can therefore be used to measure O m (e.g. White & Frenk 
1991; Fabian 1991; Briel, Henry & Bohringer 1992; White 
et al. 1993, David, Jones & Forman 1995; White & Fabian 
1995; Evrard 1997; Mohr, Mathiesen & Evrard 1999; Ettori 
& Fabian 1999; Roussel, Sadat & Blanchard 2000; Grego et 
al. 2001; Allen, Schmidt & Fabian 2002a; Allen et al. 2003, 
2004; Ettori, Tozzi & Rosati 2003; Sanderson & Ponman 
2003; Lin, Mohr & Stanford 2003; LaRoque et al. 2006). 
This method currently provides one of our best constraints 
on f2 m and is remarkably simple and robust in terms of its 
underlying assumptions. 

Measurements of the apparent evolution of the cluster 
X-ray gas mass fraction, hereafter / gas , can also be used to 
probe the acceleration of the Universe (Allen et al. 2004; see 
also Sasaki 1996, Pen 1997, Allen et al. 2002a, 2003, Ettori 
et al. 2003, LaRoque et al. 2006). This constraint originates 
from the dependence of the / gas measurements, which derive 
from the observed X-ray gas temperature and density pro- 
files, on the assumed distances to the clusters, / gas oc <f 1,5 Q 
The expectation from non-radiative hydrodynamical simu- 
lations is that for the largest (kT > 5keV), dynamically 
relaxed clusters and for measurement radii beyond the in- 
nermost core (r > J*25oo), /gas should be approximately con- 
stant with redshift (Eke et al. 1998; Crain et al. 2007). How- 
ever, possible systematic variation of / gas with redshift can 
be accounted for in a straightforward manner, so long as 
the allowed range of such variation is constrained by numer- 
ical simulations or other, complementary data (Eke et al. 
1998; Bialek, Evrard & Mohr 2001; Muanwong et al. 2002; 
Borgani et al. 2004; Kay S. et al. 2004; Ettori et al. 2004, 
2006; Kravtsov, Nagai & Vikhlinin 2005; Nagai, Vikhlinin 
& Kravtsov 2007a). 

The first clear detection of cosmic acceleration using 
the / gas technique was made by Allen et al. (2004) using 
Chandra observations of 26 hot (kT > 5keV), X-ray lumi- 
nous (Lboi ^ 10 45 /if 2 ergs -1 ), dynamically relaxed clusters 

1 To understand the origin of the / gas oc d 1 ' 5 dependence, con- 
sider a spherical region of observed angular radius 8 within which 
the mean gas mass fraction is measured. The physical size. R, is 
related to the angle 8 as R = 8d\. The X-ray luminosity emitted 
from within this region, Lxj is related to the detected flux, Fx, 
as Lx = iird^Fx, where d^ is the luminosity distance and d^ = 
o!l/(1+z) 2 is the angular diameter distance. Since the X-ray emis- 
sion is primarily due to collisional processes (bremsstrahlung and 
line emission) and is optically thin, we may also write Lx oc n 2 V, 
where n is the mean number density of colliding gas particles and 
V is the volume of the emitting region, with V = 4n(8dj\) 3 /3. 
Considering the cosmological distance dependences, we see that 
n oc d-^/d^ 5 , and that the observed gas mass within the mea- 
surement radius M gas oc nV oc d^d^ 5 . The total mass, Mtot, 
determined from the X-ray data under the assumption of hydro- 
static equilibrium, Aftot oc d\. Thus, the X-ray gas mass fraction 
measured within angle 8 is / gas = M gas /Mtot oc d^d^ 5 . 



spanning the redshift range 0.07—0.9. The total Chandra ex- 
posure used in that work, after all screening procedures were 
applied, was ~ 830ks. That study led to a ~ 3<r detection 
of the acceleration of the Universe and a tight constraint on 
the mean mass density Q m = 0.25 ± 0.04 (see also Allen et 
al. 2002a, 2003; Ettori et al. 2003; LaRoque et al. 2006) in 
excellent agreement with independent findings from CMB 
studies (e.g. Spergel et al. 2003, 2007), Type la supernovae 
(SNIa) data (e.g. Riess et al. 2004; Astier et al. 2006), galaxy 
redshift surveys (e.g. Cole et al. 2005; Eisenstein et al. 2005; 
Percival et al. 2007) and X-ray cluster number counts (e.g. 
Mantz et al. 2007). 

Here we present a significant extension of the Allen et 
al. (2004) work. Our expanded sample contains 42 clusters 
spanning the redshift range 0.05 < z < 1.1. We incorpo- 
rate new, deeper exposures for some of the original clus- 
ters, as well as new targets, approximately doubling the to- 
tal exposure time used. Our analysis method incorporates 
conservative allowances for systematic uncertainties associ- 
ated with instrument calibration, cluster physics and data 
modelling. As before, we employ rigorous selection criteria, 
restricting the analysis to the hottest, most dynamically re- 
laxed clusters. We show that this leads to remarkably small 
intrinsic scatter in the / gas measurements, with no apparent 
systematic dependence of / gas on temperature for clusters 
with kT > 5keV. Our method imposes a minimum of prior 
constraints and does not require that the density and tem- 
perature profiles of the X-ray emitting gas follow simple pa- 
rameterized forms. We make our / gas measurements for each 
cluster at the radius r2soo in the reference ACDM cosmol- 
ogy, corresponding to an angle O^soo M , f° r which the mean 
enclosed mass density is 2500 times the critical density of 
the Universe at the redshift of the cluster. This corresponds 
to about one quarter of the virial radiufl and represents a 
near-optimal choice for Chandra studies, being sufficiently 
large to provide small systematic scatter but not so large 
as to be hampered by systematic uncertainties in the back- 
ground modelling. We compare our / gas measurements to 
results from other, independent studies and to the predic- 
tions from current hydrodynamical simulations. 

Our analysis of cosmological parameters employs a 
Markov Chain Monte Carlo approach, which is efficient and 
allows for the simple inclusion of priors and a comprehensive 
study of the effects of systematic uncertainties. We present 
results based on studies of the / gas data alone (adopting 
simple priors on Qbh 2 and h) and for the / gas data used in 
combination with current CMB constraints (in which case 
the priors on Q,^h 2 and h can be dropped) and SNIa data 
(Astier et al. 2006; Riess et al. 2007; Wood-Vasey et al. 2007; 
Jha, Riess & Kirshner 2007). We highlight the power of 
the data combinations for cosmological work, particularly 
in constraining the mean matter and dark energy densities 
of the Universe and the dark energy equation of state. 

The / gas measurements are quoted for a flat ACDM 
reference cosmology with h = J/o/lOOkms -1 Mpc -1 =0.7 
and Q,m — 0.3. 



2 The virial radius is defined as the radius within which the den- 
sity contrast A c = 178 n m (z) 0,45 , with respect to the critical 
density (Lahav et al. 1991; Eke et al. 1998). 
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Table 1. Summary of the Chandra observations. Columns list the target name, observation date, detector used, observation mode, net 
exposure after all cleaning and screening processes were applied and Right Ascension (R.A.) and Declination (Dec.) for the X-ray centres. 
Where multiple observations of a single cluster have been used, these are listed separately. 
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Abell 1795(1) 


2002 Jun 10 
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oz . 




26 35 38 


Abell 1795(2) 
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ACIS-S 
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Abell 1795(3) 
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ACIS-S 
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ACIS-I 
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ACIS-S 
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2004 Sep 20 


ACIS-I 


VFAINT 


8.5 










« 


Abell 383(1) 


2000 Nov 16 


ACIS-S 


FAINT 


18.0 


02 


lo 


Uo. 


o 
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RXJ 1504. 1-0248 
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VFAINT 
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1 
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2003 Nov 18 


ACIS-S 


VFAINT 
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MACSJ0242. 6-2132 


2002 Feb 07 


ACIS-I 


VFAINT 


10.2 


02 


42 
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9 


-21 32 26 


MACSJ1427.6-2521 


2002 Jun 29 


ACIS-I 


VFAINT 


14.7 


14 


27 


39. 


4 


-25 21 02 


MACSJ2229. 8-2756 


2002 Nov 13 


ACIS-I 


VFAINT 


11.8 


22 


29 


45. 


3 


-27 55 37 


MACS J0947. 2+7623 


2000 Oct 20 


ACIS-I 


VFAINT 


9.6 


09 


47 
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1 


76 23 14 


MACSJ1931. 8-2635 


2002 Oct 20 


ACIS-I 


VFAINT 


12.2 


19 


31 
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6 
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MACSJ1115.8+0129 
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VFAINT 
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11 


15 
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1 
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2 X-RAY OBSERVATIONS AND ANALYSIS 
2.1 Sample selection 

Our sample consists of 42 hot, X-ray luminous, dynami- 
cally relaxed galaxy clusters spanning the redshift range 
0.05 < z < 1.1. The systems have mass weighted X-ray 
temperatures measured within r2soo, fcTbsoo 5keV and 
exhibit a high degree of dynamical relaxation in their Chan- 
dra images (Million et al. 2007, in prep.), with sharp cen- 
tral X-ray surface brightness peaks, short central cooling 
times (t coo i < a few 10 9 yr) minimal isophote centroid vari- 
ations (e.g. Mohr et al. 1995) and low X-ray power ra- 
tios (Buote & Tsai 1995, 1996; Jeltema et al. 2005). Al- 
though target selection is based only on these morphological 
X-ray characteristics, the clusters also exhibit other signa- 
tures of dynamical relaxation including minimal evidence 
for departures from hydrostatic equilibrium in X-ray pres- 
sure maps (Million et al. 2007, in prep.). The notable excep- 
tions are Abell 2390, RXJ1347.5-1145, MACS1427.3+4408 
and MACSJ0744.9+3927, for which clear substructure is ob- 
served between position angles of 255-15 degrees, 90-190 
degrees, 160-280 degrees and 210-330 degrees, respectively 



(Allen, Schmidt & Fabian 2002b; Morris et al. 2007, in prep.; 
Million et al. 2007, in prep.). The regions associated with ob- 
vious substructure in these clusters have been excluded from 
the analysis. The bulk of the clusters at z > 0.3 were iden- 
tified in the MACS survey (Ebeling, Edge & Henry 2001; 
Ebeling et al. 2007). Of the 70 MACS clusters with suffi- 
cient data on the Chandra archive at the time of writing 
to enable detailed spatially-resolved spectroscopy, 22/70 are 
identified as being sufficiently relaxed to be included in the 
present study. 

The restriction to clusters with the highest possible de- 
gree of dynamical relaxation, for which the assumption of 
hydrostatic equilibrium should be most valid, minimizes sys- 
tematic scatter in the / gas data (Section 15. 3[) and allows for 
the most precise and robust determination of cosmological 
parameters. The restriction to the hottest (kT > 5keV), re- 
laxed systems further simplifies the analysis: for galaxies, 
groups and clusters with kT < 4keV, the baryonic mass 
fraction is both expected and observed to rise systemati- 
cally with increasing temperature, with the systematic scat- 
ter being largest in the coolest systems (e.g. Bialek et al. 
2001; Muanwong et al. 2002; Ettori et al. 2004; Kravtsov, 
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Table 1. Summary of the Chandra observations - continued 
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12 
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8.70 


11 


27 
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.2 


11 


07 
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2005 


Apr 
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MACSJ1423.8+2404 
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58. 
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ACIS-I 


VFAINT 
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11. 
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12 
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15.4 


07 


11 


17. 


.5 


37 
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17 



et al. 2005; Vikhlinin et al. 2006). As shown in Sections ED 
and 15.31 for the hot, relaxed clusters studied here, f gas ex- 
hibits no dependence on temperature and the intrinsic scat- 
ter is small. 



2.2 Data reduction 

The Chandra observations were carried out using the Ad- 
vanced CCD Imaging Spectrometer (ACIS) between 1999 
August 30 and 2005 June 28. The standard level-1 event 
lists produced by the Chandra pipeline processing were re- 
processed using the CI AO (version 3.2.2) software package, 
including the appropriate gain maps and calibration prod- 
ucts. Bad pixels were removed and standard grade selections 
applied. Where possible, the extra information available in 
VFAINT mode was used to improve the rejection of cos- 
mic ray events. The data were cleaned to remove periods 
of anomalously high background using the standard energy 
ranges and time bins recommended by the Chandra X-ray 
Center. The net exposure times after cleaning are summa- 
rized in Table [T] The total good exposure is 1.63 Ms, ap- 
proximately twice that of the Allen et al. (2004) study. 



2.3 Spectral analysis 

The spectral analysis was carried out using an updated ver- 
sion of the techniques described by Allen et al. (2004) and 
Schmidt & Allen (2007). In brief, concentric annular spec- 
tra were extracted from the cleaned event lists, centred on 
the coordinates listed in Table [T] Emission associated with 
X-ray point sources or obvious substructure (Table [2} was 
excluded. The spectra were analysed using XSPEC (ver- 
sion 11.3; Arnaud 1996), the MEKAL plasma emission code 
(Kaastra & Mewe 1993; incorporating the Fe-L calculations 
of Liedhal, Osterheld & Goldstein 1995) and the photoelec- 
tric absorption models of Balucinska-Church & McCammon 
(1992). The emission from each spherical shell was modelled 
as a single phase plasma. The abundances of the elements 
in each shell were assumed to vary with a common ratio, Z, 
with respect to Solar values. The absorbing column densities 
were fixed to the Galactic values determined from HI studies 
(Dickey & Lockman 1990), with the exception of Abell 478 
and PKS0745-191 where the value was allowed to fit freely. 
(For Abell 478, the absorbing column density was allowed 
to vary as a function of radius, as was shown to be required 
by Allen et al. 1993). We have included standard correc- 
tion factors to account for time-dependent contamination 
along the instrument light path. In addition, we have incor- 
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porated a small correction to the High Resolution Mirror 
Assembly model in CIAO 3.2.2, which takes the form of an 
'inverse' edge with an energy, E=2.08 keV and optical depth 
r — —0.1 (H. Marshall, private communication) and also 
boosted the overall effective area by six per cent, to better 
match later calibration data (A. Vikhlinin, private commu- 
nication) . These corrections lead to an excellent match with 
results based on later calibration data, available in CIAO 
3.4. Only data in the 0.8 — 7.0 keV energy range were used 
in the analysis (with the exceptions of the earliest observa- 
tions of 3C 295, Abell 1835 and Abell 2029, where a wider 
0.6 to 7.0 keV band was used to enable better modelling of 
the soft X-ray background). 

For the nearer clusters (z < 0.3), background spectra 
were extracted from the blank-field data sets available from 
the Chandra X-ray Center. These were cleaned in an iden- 
tical manner to the target observations. In each case, the 
normalizations of the background files were scaled to match 
the count rates in the target observations measured in the 
9.5-12keV band. Where required, e.g. due to the presence 
of strong excess soft emission in the field, a spectral model 
for additional soft background emission was included in the 
analysis. For the more distant systems (as well as for the first 
observation of Abell 1835, the ACIS-I observation of Abell 
383, and the observations of Abell 2537, RXJ 2129.6+0005 
and Zwicky 3146) background spectra were extracted from 
appropriate, source free regions of the target data sets. (We 
have confirmed that similar results are obtained using the 
blank-field background data sets.) In order to minimize sys- 
tematic errors, we have restricted our spectral analysis to 
radii within which systematic uncertainties in the back- 
ground subtraction (established by the comparison of differ- 
ent background subtraction methods) are smaller than the 
statistical uncertainties in the results. All results are drawn 
from ACIS chips 0,1,2,3 and 7 which have the most accurate 
calibration, although ACIS chip 5 was also used to study the 
soft X-ray background in ACIS-S observations. 

Separate photon-weighted response matrices and effec- 
tive area files were constructed for each region using cali- 
bration files appropriate for the period of observations. The 
spectra for all annuli for a given cluster were modelled simul- 
taneously in order to determine the deprojected X-ray gas 
temperature and metallicity profiles, under the assumption 
of spherical symmetry. The extended C-statistic, available 
in XSPEC, was used for all spectral fitting. 

2.4 Measuring the mass profiles 

The details of the mass analysis and results on the total mass 
and dark matter profiles are presented by Schmidt & Allen 
(2007). In brief, X-ray surface brightness profiles in the 0.8- 
7.0keV band were extracted from background subtracted, 
flat-fielded Chandra images with 0.984 x 0.984arcsec 2 pix- 
els. The profiles were centered on the coordinates listed in 
Table [T] Under the assumptions of hydrostatic equilibrium 
and spherical symmetry, the observed X-ray surface bright- 
ness profiles and deprojected X-ray gas temperature profiles 
may together be used to determine the X-ray emitting gas 
mass and total mass profiles in the clusters. For this anal- 
ysis, we have used an enhanced version of the Cambridge 
X-ray deprojection code described by e.g. White, Jones & 
Forman (1997). This method is particularly well suited to 



the present task in that it does not use parametric fitting 
functions for the X-ray temperature, gas density or surface 
brightness in measuring the mass; the use of such functions 
introduces strong priors that complicate the interpretation 
of results and, in particular, can lead to an underestimation 
of uncertainties. The only additional important assumption 
in the analysis is the choice of a Navarro, Frenk & White 
(1995, 1997; hereafter NFW) model to parameterize the to- 
tal (luminous-plus-dark) mass distributions: 



p(r) = 



p c (z)S c 



(r/r s ) (l + r/r B y 



(1) 



where p(r) is the mass density, p c (z) = 3H(z) 2 /8ivG is the 
critical density for closure at redshift z, r B is the scale ra- 
dius, c is the concentration parameter (with c = r2oo/Vs) 
and S c = 200c 3 /3 [ln(l + c) - c/(l + c)]0 Schmidt & Allen 
(2007) show that the NFW model provides a good descrip- 
tion of the mass distributions in the clusters studied here. 

Given the observed surface brightness profile and a par- 
ticular choice of parameters for the total mass profile, the 
deprojection code is used to predict the temperature profile 
of the X-ray gas. (In detail, the median model tempera- 
ture profile determined from 100 Monte-Carlo simulations 
for each mass model is used.) This model temperature pro- 
file is then compared with the observed spectral, deprojected 
temperature profile and the goodness of fit is calculated us- 
ing the sum over all temperature bins 



* 2 = E 

all bins 



'Tob 



(5 obs 



(2) 



where T bs is the observed, spectral deprojected tempera- 
ture profile and T mo dd is the model, rebinned to the same 
spatial scale. For each cluster, the mass parameters are 
stepped over a grid of values and the best-fit values and un- 
certainties determined via % 2 minimization techniques. The 
X-ray emitting gas density, pressure, entropy, cooling time 
and mass, and the integrated X-ray gas mass fraction, / gas , 
are then determined in a straightforward manner from the 
Monte-Carlo simulations and \ 2 values at each grid point. 

A number of systematic issues affect the accuracy of 
the /g as measurements and their interpretation; these are 
discussed in detail in Section ^. 21 In particular, our analysis 
incorporates allowances for effects associated with calibra- 
tion and modelling uncertainties and non-thermal pressure 
support in the X-ray emitting gas, employing priors that 
span conservative ranges for the likely magnitudes of these 
effects. 

Finally, for a number of the clusters, noticeable sub- 
structure is present at small radii. This is likely to result 
from interactions between the central radio sources and sur- 
rounding gas {e.g. Bohringer et al. 1993; Fabian et al. 2000, 
2003a, 2005, 2006; Birzan et al. 2004; Dunn & Fabian 2004; 
Forman et al. 2005; Dunn, Fabian & Taylor 2005; Allen et 



3 Note that the outermost pressure, at the limit of the X-ray 
surface brightness profile, is fixed using an iterative method that 
ensures a smooth, power law pressure gradient in these regions. 
The model temperature profiles, for radii spanned by the spectral 
data, are not sensitive to any reasonable choices for the outer 
pressures. 
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Table 2. Clusters with regions of localized substructure that 
have been excluded or down- weighted in the analysis. Column 
two lists the position angles (PA) that have been excluded in 
the case of Abell 2390, RXJ1347.5-1145, MACS1427.3+4408 and 
MACSJ0744.9+3927. Column 3 lists the radii (in /i^kpc) within 
which the spectral data have been down-weighted by including a 
systematic uncertainty of ±30 per cent in quadrature with the 
statistical errors on the temperature measurements. 



Cluster Excluded P.A. Down-weighted r 



Abell 1795 75 

Abell 2029 30 

Abell 478 15 

PKS0745-191 55 

Abell 1413 40 

Abell 2204 75 

Abell 383 40 

RXJ1504. 1-0248 80 

Abell 2390 255 - 15 50 

RXJ2129.6+0005 40 

Zwicky 3146 240 

Abell 2537 40 

MACSJ2229.8-2756 40 

MACSJ0947.2+7623 40 

MACSJ1931. 8-2635 40 

MACSJ1115.8+0129 85 

MACSJ1532.9+3021 40 

RXJ1347.5-1145 90 - 190 

MACSJ1621. 6+3810 45 

MACSJ1427.3+4408 160 - 280 

MACSJ0744.9+3927 210 - 330 



al. 2006; Rafferty et al. 2006) and/or 'sloshing' of the X-ray 
emitting gas within the central potentials (e.g. Churazov et 
al. 2003; Markevitch et al. 2003; Ascasibar & Markevitch 
2006). The regions affected by such substructure are listed 
in Table [5] A systematic uncertainty of ±30 per cent has 
been added in quadrature to all spectral results determined 
from these regions, leading to them having little weight in 
the mass analysis. 

2.5 The stellar baryonic mass fraction 

Observations of nearby and intermediate redshift clusters 
show that for clusters in the mass/temperature range stud- 
ied here, the average mass fraction in stars (in galaxies 
and intracluster light combined) / star ~ 0.16ft7o 5 / gas (Lin 
& Mohr 2004; see also White et al. 1993; Fukugita, Hogan 
& Peebles 1998; Balogh et al. 2001). 

For the present analysis, we ideally require the ratio 
s — /star //gas measured within r25oo for each cluster. How- 
ever, such measurements are not yet available for the bulk 
of the clusters studied here. For hot, massive clusters, the 
relative contribution of the central dominant galaxy to the 
overall cluster light is less significant than for cooler, less 
massive systems {e.g. Lin & Mohr 2004). We have therefore 
assumed that the stellar mass fraction within r25oo is similar 
to that measured within the virial radius i.e. s = O.I6/175 5 , 
but have both included a conservative 30 per cent Gaus- 
sian uncertainty in this value and allowed for evolution at 
the ±20 per cent level, per unit redshift interval. Since the 



stellar mass accounts for only ~ 14 per cent of the over- 
all baryon budget within r2soo and less than 2 per cent of 
the total mass, these systematic uncertainties do not have a 
large effect on the overall error budget. A program to mea- 
sure the evolution of the optical baryonic mass content of 
the largest relaxed clusters is underway. 



3 THE X-RAY GAS MASS FRACTION 
MEASUREMENTS 

3.1 New /g as measurements 

As mentioned above, in compiling the results on the X-ray 
gas mass fraction, / gas , we have adopted a canonical mea- 
surement radius of r25oo- The r2soo value for each cluster is 
determined directly from the Chandra data, with confidence 
limits calculated from the \ 2 grids. In general, the values 
are well-matched to the outermost radii at which reliable 
temperature measurements can be made from the Chan- 
dra data, given systematic uncertainties associated with the 
background modelling. 

Fig. [TJa) shows the observed /gas(r) profiles for the 
six lowest redshift clusters in the sample, for the reference 
ACDM cosmology. Although some dispersion in the profiles 
is present, particularly at small radii, the profiles tend to- 
wards a common value at r25oo- Fitting the / gas measure- 
ments at r25oo for the six lowest-redshift systems with a 
constant value we obtain / gas = 0.113±0.003, with \ 2 = 4.3 
for 5 degrees of freedom. Fitting the results for all 42 clusters 
gives / gas = 0.1104 ± 0.0016, with \ 2 = 43.5 for 41 degrees 
of freedom. 

Fig- Hlb) shows the / g as(»"/?"250o) profiles for all 42 clus- 
ters in the sample. Fitting the data in the range 0.7 — 
1.2r25oo with a power-law model, we measure /g as = 0.1105± 
0.0005(r/r250o)°- 214±0 022 . Note that the error bars on the 
mean / gas measurements quoted above reflect only the sta- 
tistical uncertainties in these values. A systematic uncer- 
tainty of ~ 10 — 15 per cent in the global, absolute / gas nor- 
malization is also present due to uncertainties in e.g. instru- 
ment calibration, X-ray modelling and non-thermal pressure 
support; this must be accounted for in the determination of 
cosmological constraints (Section 4.2). 

Table [3] summarizes the results on the X-ray gas mass 
fraction for each cluster measured at r2soo, together with 
the T2500 values, for the reference ACDM cosmology. Fig. 
[5] shows a comparison of the / gas results, plotted as a func- 
tion of redshift, for the reference ACDM cosmology and a 
flat, standard cold dark matter (SCDM) cosmology with 
Sl m = 1.0, h = 0.5. Whereas the results for the ACDM 
cosmology appear consistent with the expectation of a con- 
stant /g a s (z) value from non-radiative simulations (e.g. Eke 
et al. 1998; Crain et al. 2007), as evidenced by the accept- 
able \ 2 value quoted above, the results for the reference 
SCDM cosmology indicate a clear, apparent drop in / gas 
as the redshift increases. The \ 2 value obtained from a fit 
to the SCDM data with a constant model, \ 2 = 144 for 
41 degrees of freedom, shows that the SCDM cosmology is 
clearly inconsistent with a prediction that f gas (z) should be 
constant. 

Table [3] also lists the mass- weighted temperatures mea- 
sured within r25oo for each cluster. Fig. [3] shows / gas as a 
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Figure 1. The X-ray gas mass fraction profiles for the ACDM reference cosmology (C m = 0.3, Qa = 0.7, h = 0.7) with the radial axes 
scaled in units of r25oo- ( a: left panel) Results for the six lowest redshift clusters with z < 0.15 (b: right panel) Results for the entire 
sample. Note f ga , B (r) is an integrated quantity and so error bars on neighbouring points in a profile are correlated. 



o 








0.5 




Figure 2. The apparent variation of the X-ray gas mass fraction measured within r2500 as a function of redshift for the (a: left panel) 
reference ACDM and (b: right panel) reference SCDM (f2 m = 1.0, Oa = 0.0, h = 0.5) cosmologies. The plotted error bars are statistical 
root-mean-squarc lcr uncertainties. The global, absolute normalization of the / gas values should be regarded as uncertain at the ~ 10 — 15 
per cent level due to systematic uncertainties in instrument calibration, modelling and the level of non-thermal pressure support (Section 
4.2). 



function of fcTbsoo for the reference ACDM cosmology. The 
dotted line in the figure shows the best-fitting power law 
model, /gas(V25oo) x fc^aloO) which provides a good descrip- 
tion of the data (\ 2 = 43.5 for 40 degrees of freedom) and 
is consistent with a constant value (a — 0.005 ± 0.058). The 
solid lines show the 2a limits on the steepest and shallowest 
allowed power law models. It is clear from the figure that 
/gas is independent of temperature for the clusters in the 
present sample. 

3.2 Comparison with previous / gas results 

Approximately 0.75 Ms of the ~ 1.6 Ms of Chandra data 
used here were also included in the Allen et al. (2004) study. 
The current work includes a re-analysis of those data us- 



ing improved calibration information, where available. The 
/ gaa results from the two studies show excellent overall agree- 
ment: the new / ga s values are, on average, ~ 6 per cent 
lower than those reported by Allen et al. (2004), a differ- 
ence consistent with expectations given the modification to 
the effective area calibration described in Section \2. 31 

LaRoque et al. (2006) present / gas measurements for 38 
X-ray luminous clusters, including 10 of the large, dynami- 
cally relaxed systems studied here. Their best-fit results at 
»"2500 are in good overall agreement with the present work, 
with their / gas values being, on average, ~ 6 per cent higher 
than those reported here, for the systems in common. 

Pointecouteau et al. (2004) present an analysis of 
XMM-Newton data for Abell 478, for which they measure 
an / g as value at r25oo of 0.13 ±0.02, in good agreement with 
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Table 3. The redshifts, r25oo values, mean mass-weighted temperatures within r2500i an d the X-ray gas mass fractions within r2soo for 
the reference ACDM cosmology. Error bars are statistical uncertainties and are quoted at the 68 per cent confidence level. A systematic 
uncertainty of ~ 10 — 15 per cent is associated with the global, absolute normalization of the / gas values due to uncertainties in instrument 
calibration, X-ray modelling and non-thermal pressure support (Section 4.2). The redshifts for the MACS clusters are from Ebeling et 
al. (2007, in preparation). 
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this work. These authors also report a value of 0.11 for Abell 
1413, based on the data of Pratt & Arnaud (2002), which is 
consistent with the results reported here. 

Vikhlinin et al. (2006) present / gas measurements for 
13 clusters of which six are in common with this study. On 
average, the Vikhlinin et al. (2006) / gas results are ~ 10 per 
cent lower than those reported here after correcting their 
values to the same reference ACDM cosmology. 

We note that the statistical uncertainties on the 



/ gas measurements listed in Table [3] are, typically, larger 
than those reported by other authors. Two contributing fac- 
tors to this difference are: 1) that the present analysis does 
not impose strong priors on the shapes of the temperature 
and density profiles in the clusters through the use of para- 
metric models (the use of such parameterizations can lead 
to spuriously tight constraints in cases where they do not 
provide an adequate description of the data); and 2) the 
/ gaa measurement errors reported here are marginalized over 
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Figure 3. The X-ray gas mass fraction as a function of mass- 
weighted temperature measured within r2500 for the reference 
ACDM cosmology. The dotted line shows the best-fitting power 
law model which provides a good description of the data (v 2 = 
43.5 for 40 degrees of freedom) and is consistent with a constant 
value (slope a = 0.005±0.058). The solid lines show the 2<r limits 
on the slopes allowed by the data. The figure demonstrates that 
/gas is essentially independent of temperature for the massive, 
dynamically relaxed clusters in the present sample. 



the uncertainties in all other parameters, including the un- 
certainties in r 2 5oo. 



4 COSMOLOGICAL ANALYSIS 

4.1 Markov Chain Monte Carlo method 

Our determination of cosmological parameters uses a 
Markov Chain Monte Carlo method. We employ a modified 
version of the CosmoMC cod^fl of Lewis & Bridle (2002; see 
Rapetti et al. 2005, 2007 for details of the enhancements), 
which uses a Metropolis-Hastings Markov Chain Monte 
Carlo (MCMC) algorithm to explore parameter space. We 
run the code on four to sixteen processors simultaneously, 
creating multiple chains and using the Message Passing In- 
terface (MPI) to dynamically update the proposal matrix 
based on the covariance of post-burn-in samples. This leads 
to a much faster convergence than would be obtained from 
a single chain run on a single compute node. 

Convergence is assessed using the Gelman-Rubin cri- 
terion (Gelman & Rubin 1992). Convergence is deemed ac- 
ceptable when the ratio of between-chain to mean-chain vari- 
ances, R, satisfies R — 1 < 0.1. (We have also visually com- 
pared individual chains to ensure that consistent final results 
were obtained.) In general, our combined chains typically 
have lengths of at least 10 5 samples and have R — 1 << 0.1. 
(For the evolving-io models, R— 1 ~ 0.1.) Conservative burn- 
in periods of at least 10000 samples were allowed for each 
chain. 



http: / /cosmologist.info/cosmomc/ 



4.2 Analysis of the / gas data: modelling and 
systematic allowances 

The differences between the shapes of the / gas (z) curves in 
Figs. [21 a) and (b) reflect the dependence of the measured 
/gas values on the assumed angular diameter distances to the 
clusters. Under the assumption (Section 1) that / gas should, 
in reality, be approximately constant with redshift, as sug- 
gested by non-radiative simulations of large clusters (Eke et 
al. 1998; Crain et al. 2007; uncertainties in the predictions 
from simulations are discussed below) inspection of Fig. [2] 
would clearly favour the ACDM over the SCDM cosmology. 

To determine constraints on cosmological parameters, 
it is not necessary to generate f sas (z) data sets for every 
cosmology of interest and compare them to the expected 
behaviour. Rather, one can fit a single, reference f sas (z) data 
set with a model that accounts for the expected apparent 
variation in / gas (z) as the underlying cosmology is varied. 
We choose to work with the ACDM reference cosmology, 
although similar results can in principle be derived for other 
reference cosmologies. 

The model fitted to the reference ACDM data is 



KA-yb(z) (Q. 
l + s(z) 



/Ob 

Vftm 



4 CDM ( 2 ) 



d\(z) 



(3) 



where dx{z) and d^ CDM (z) are the angular diameter dis- 
tances to the clusters in the current test model and reference 
cosmologies, 



d < 



H (l + z)V£h 



sinh 



dz 



(4) 



with E{z) defined as in Section T4.4I The factor A in Equa- 
tion [3] accounts for the change in angle subtended by r2soo 
as the underlying cosmology is varied 



A = 



gACDM 
'2500 



025 



H(z)d A (z) 
[H(z)d A (z)Y 



(•») 



Here, r\ is the slope of the f gas (r /r25oo) data in the region of 
^2500, as measured for the reference ACDM cosmology. For 
simplicity, we use the best-fit average slope of 77 = 0.214 ± 
0.022 determined from a fit to the whole sample over the 
range 0.7 < r/r 2 5oo < 1.2 (Section O and marginalize over 
the slope uncertainty. This angular correction factor, which 
is close to unity for all cosmologies and redshifts of interest, 
has not been employed in previous studies and, indeed, can 

5 To see the origin of the correction factor A, recall that equa- 
tion[3]predicts the / gas value at the measurement radius in the ref- 
erence ACDM cosmology. This measurement radius corresponds 
to a fixed angle 025O? M f° r eacn cluster, which will differ slightly 
from 2 5oo, the angle corresponding to r25oo for that cluster in 
the current test cosmology. The mass contained within radius 
7*2500, M 25 oo = Iu 4 7rr3 500 p c r it/3. Given that the temperature, 
and temperature and density gradients, in the region of 6*2500 
are likely to be approximately constant, the hydrostatic equation 
gives M2500 25 r2500- Thus, since p cr it = 3H(z) 2 /8nG, we have 
r 2500 25 H(z) — 1 , and the angle spanned by r25oo a t redshift z, 
6*2500 = f250o/dA 25 (H(z)d A )~ 1 . Since the / gas profiles follow 
a smooth power law form in the region of 6*2500, the ratio of the 
model / gas value at O^SS^ *° that at 6*2500 can be described by 
Equation [5] 
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be neglected without significant loss of accuracy for most 
work. Nevertheless, we include it here for completeness and 
note that its inclusion leads to slightly tighter constraints 
on dark energy than would otherwise be obtained. 

The parameter 7 in equation [3] models non-thermal 
pressure support in the clusters. Based on hydrodynami- 
cal simulations, Nagai et al. (2007a) estimate a bias of ~ 9 
per cent in / gas measurements at r2soo for relaxed clusters. 
This bias originates primarily from subsonic motions in the 
intracluster gas and, as discussed by those authors (see also 
Section |5.3|I . can be regarded as an upper limit, given ob- 
servational indications that the gas viscosity in real clusters 
appears likely to exceed that modelled in the simulations. 
For the large, relaxed clusters and measurement radii of 
interest here, non-thermal pressure support due to cosmic 
rays (Pfrommer et al. 2007) and magnetic fields (Dolag & 
Schindler 2000) is expected to be small. Based on these con- 
siderations, our default analysis assumes a uniform prior of 
1.0 < 7 < 1.1, although we also consider the case where the 
non-thermal pressure support may be up to twice as large 
i.e. 1.0 < 7 < 1.2. 

The parameter s(z) = sq(1 + ot B z) in equation [3] models 
the baryonic mass fraction in stars. As discussed in Sec- 
tion 12.51 we include a 30 per cent Gaussian uncertainty on 
so, such that so = (0.16 ± 0.05)/i7 5 i an< A a 20 per cent uni- 
form prior on ct s , such that —0.2 < a s < 0.2, allowing for 
evolution in the stellar baryonic mass fraction of ±20 per 
cent per unit redshift interval. 

The factor b(z) = bo(l + a b z) is the 'depletion' or 'bias' 
factor i.e. the ratio by which the baryon fraction measured 
at T2500 is depleted with respect to the universal mean; such 
depletion is a natural consequence of the thermodynamic 
history of the gas. The non-radiative simulations of hot, 
massive clusters published by Eke et al. (1998; see also Crain 
et al. 2007) give bo — 0.83 ± 0.04 at r2soo, and are consis- 
tent with no redshift evolution in 6 for z < 1. We use these 
simulations as a benchmark because other simulations that 
include cooling currently tend to significantly over-produce 
young stars in the largest galaxies (see e.g. Balogh et al. 
2001), which is problematic for the prediction of b(z). We 
note also the good agreement between the observed, scaled 
/gas(r) profiles determined from the Chandra data and the 
b(r) profiles for the three most relaxed clusters in the sim- 
ulations of Eke et al. (1998; the red curves in Fig [J}; this 
suggests that the non-radiative simulations provide a use- 
ful approximation for the purpose of predicting b(z). (The 
profiles for the less relaxed simulated clusters are shown as 
dashed green curves in the figure.) Nevertheless, to account 
for systematic uncertainties in the predictions of b(z), we 
include a conservative 20 per cent uniform prior on bo, such 
that 0.65 < 60 < 1-0, and allow for moderate, systematic 
evolution in b(z) over the observed redshift range, setting 
—0.1 < at < 0.1. This encompasses a range of evolution 
allowed by recent simulations including various approxima- 
tions to the detailed baryonic physics {e.g. Kay et al. 2004, 
Ettori et al. 2006, Crain et al. 2007, Nagai et al. 2007a). 

The factor K in equation [3] is a 'calibration' constant 
that parameterizes residual uncertainty in the accuracy of 
the instrument calibration and X-ray modelling. Contribut- 
ing factors include uncertainty in the instrument effective 
area, variations in element abundance ratios, modelling the 
effects of gas clumping and asphericity (the latter effects are 



expected to be small for large, relaxed clusters; Nagai et al. 
2007a. See also Piffaretti, Jetzer & Schindler 2003, Gavazzi 
2005). We conservatively include a 10 per cent Gaussian un- 
certainty in K to model the combined effect of these factors, 
such that K = 1.0 ± 0.1. The small intrinsic dispersion in 
/ gas values (Section I5.3[l means that Malmquist bias is ex- 
pected to have a negligible effect on the derived cosmological 
parameters. Uncertainties associated with other systematic 
factors are expected to be negligible in comparison to the 
allowances listed above. 

In cases where the Chandra / gas data are not combined 
with CMB data, we include simple Gaussian priors on Q b h 2 
and h. Two separate sets of priors were used: 'standard' 
priors with Q. h h 2 = 0.0214 ± 0.0020 (Kirkman et al. 2003) 
and h = 0.72 ± 0.08 (Freedman et al. 2001), and 'weak' 
priors in which the nominal uncertainties were tripled to 
give Q b h 2 = 0.0214 ± 0.0060 and h = 0.72 ± 0.24. In cases 
where the CMB data are included, no priors on Q b h 2 or h 
are needed or used. The complete set of standard priors and 
allowances included in the / gas analysis are summarized in 
Table H 

Finally, we note how inspection of equation [3] can pro- 
vide useful insight into the strength of the / gas experiment. 
The pre-factors before the square brackets shows how the 
normalization of the / gaa (z) curve is used to constrain f2 m , 
given prior information on fib, h, K, 7, b and s. The ratio of 
distances inside the square brackets (and to a small extent 
the angular correction factor) shows how the shape of the 
/ gas (z) curve constrains the geometry of the Universe and 
therefore dark energy. The combination of information from 
both the normalization and shape breaks the degeneracy be- 
tween f2 m and the dark energy parameters in the distance 
equations. 

4.3 Other data used in the analysis 

In addition to the analysis of the Chandra / gas data alone, 
we have examined the improved constraints on cosmological 
parameters that can be obtained through combination of the 
/ gas data with CMB and SNIa studies. 

Our analysis of CMB observations uses the three- 
year WMAP temperature (TT) data for multipoles 
I < 1000 (Hinshaw et al. 2007; Spergel et al. 2007) 
and temperature-polarization (TE) data for I < 450 
(Page et al. 2007). We use the October 2006 ver- 
sion of the WMAP likelihood code available from 
I http :/ /lambda.gsfc.nasa. gov /product/map / current /m_ sw. cfm\ 
Like most authors, we have ignored the small contribution to 
the TT data expected to arise from the Sunyaev-Z'eldovich 
(SZ) effect in clusters and groups (e.g. Komatsu & Seljak 

2002) and do not account for gravitational lensing of the 
CMB (Lewis & Challinor 2006), which has a negligible 
effect on the derived cosmological parameters. To extend 
the analysis to higher multipoles (smaller scales), we also 
include data from the Cosmic Background Imager (CBI; 
Mason et al. 2003; Pearson et al. 2003), the Arcminute 
Cosmology Bolometer Array Receiver (ACBAR; Kuo et al. 

2003) and BOOMERanG (Jones et al. 2006; Montroy et al. 
2006; Piacentini et al. 2005), as incorporated into the cur- 
rent version of the CosmoMC code (Lewis & Bridle 2002). 
We use a modified version of CAMB (Lewis, Challinor 
& Lasenby 2000) to calculate CMB power spectra, which 
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Figure 4. The X-ray depletion or bias factor, b (i.e. the enclosed 
baryon fraction relative to the universal value) as a function of ra- 
dius, in units of the virial radius r v ; r , from the simulations of Eke 
et al. (1998). The simulated clusters have similar masses to the 
systems studied here. The results (at zero redshift) for the three 
most dynamically relaxed clusters in the simulations are shown 
as bold red curves. Less relaxed simulated clusters are shown as 
dashed green curves. The Chandra observations for the six low- 
est redshift clusters in the / gas sample are plotted as blue cir- 
cles, with error bars. (The Chandra profiles are identical to those 
shown in Fig[T] but are scaled assuming Sl m = 0.27, Cb = 0.0413 
and r25oo = 0.25 r v ; r .) The agreement between the observed and 
predicted profiles argues that the non-radiative simulations pro- 
vide a reasonable approximation for the purpose of predicting the 
baryonic mass distributions. 

includes a consistent treatment of the effects of dark energy 
perturbations for evolving— w models (Rapetti et al. 2005; 
we assume that the sound speed in the dark energy fluid is 
equal to the speed of light). 

Our analysis of SNIa data uses two separate supernova 
samples. In the first case, we use the compilation of Davis et 
al. (2007) which includes results from the ESSENCE survey 
(60 targets; Wood-Vasey et al. 2007, Miknaitis et al. 2007), 
the SNLS first year data (57 targets; Astier et al. 2006), 45 
nearby supernovae (Jha et al. 2007) and the 30 high-redshift 
supernovae discovered by HST and reported by Riess et al. 
(2007) for which a 'gold' rating was awarded. This sample 
includes 192 SNIa in total. The second supernova sample 
is the full 'gold' sample of Riess et al. (2007) which totals 
182 SNIa, including the HST-discovered objects. For both 
samples we marginalize analytically over the absolute nor- 
malization of the distance moduli. 

4.4 Dark Energy models 

We have considered three separate dark energy models in 
the analysis: 1) standard ACDM, for which the dark en- 
ergy equation of state w = —1; 2) a model that allows any 
constant dark energy equation of state, including 'phantom' 
models with w < — 1; 3) a model in which the dark energy 
equation of state is allowed to evolve as 

w e tZ + woZt w ct (1 — a)a t + wo(l — a t )a . . 

w = = — (o) 

z + z t a(l-2a t )+a t 



where a — 1/(1 + z) is the scale factor, wo and w ct are the 
equation of state at late (present day) and early times, and 
Zt and a t are the redshift and scale factor at the transi- 
tion between the two, respectively (Rapetti et al. 2005; see 
also Chevallier & Polarski 2001; Linder 2003; Corasaniti et 
al. 2003; Linder 2007). We employ a uniform prior on the 
transition scale factor such that 0.5 < a t < 0.95. As dis- 
cussed by Rapetti et al. (2005), this model is both more 
general and more applicable to current data, which primar- 
ily constrain the properties of dark energy at redshifts z < 1, 
than models which impose a transition redshift z = 1, e.g. 
w(a) = wo + w a (l — a). 

Energy conservation of the dark energy fluid leads to 
an evolution of the energy density with scale factor 

/ \ -3-3 f " ™ {a , - da ' 

Pde(a) = pde.oa e J i a , (7) 

where pdc,o is the energy density of the dark energy fluid 
today. Using the parameterization of equation (JS| we obtain 

/ — da' — w ct In a + (w ct — w )g(a; a t ) , (8) 

Ji a 

with 

g{a . at ) = ( l n ( i-"' . ( 9 ) 

yy ' ; U-2ot/ \a(l - 2a t ) + at J 

The Friedmann equation, which relates the first time 
derivative of the scale factor of the Universe to 
the total density, can be conveniently expressed as 
{a/a) 2 = H{a) 2 = HiE(a) 2 , with 

E{a) = y/n m a- 3 + n DE /(a) + Q k a~ 2 . (10) 

Here fik is the curvature, Ode is the dark energy density 
and f(a) is its redshift dependence. (Note that we have ig- 
nored the density contributions from radiation and relativis- 
tic matter in this expression, although they are included in 
the analysis.) For our most general dark energy parameter- 
ization (Equation [6]) 

fig} — o -3(l-H"<rt) e -3(t0et— s"o)s(«;"t) (--]^ 

For ACDM cosmologies, the dark energy density is constant 
and f(a) = 1. For w < — 1 the dark energy density increases 
with time. For constant w models with w < —1/3, dark 
energy accelerates the expansion of the universe. (The re- 
sults from a purely kinematic modelling of the data, which 
does not rely on the Friedmann equation and is independent 
of the assumptions of General Relativity, are discussed by 
Rapetti et al. 2007). 

Our combined analysis of Chandra / gas , SNIa and CMB 
data therefore has up to ten interesting parameters: the 
physical dark matter and baryon densities in units of the 
critical density, the curvature Q^, the ratio of the sound 
horizon to the angular diameter distance for the CMB 
(Kosowsky, Milosavljevic & Jimenez 2002), the amplitude 
of the scalar power spectrum, the scalar spectral index, the 
optical depth to reionization, and up to three parameters 
associated with the dark energy equation of state: wo, w ct 
and at. In all cases, we assume an absence of both tensor 
components and massive neutrinos and, for the analysis of 
the CMB data alone, include a wide uniform prior on the 
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Table 4. Summary of the standard systematic allowances and 
priors included in the Chandra / gas analysis. The priors on Q,hh 2 
and h (Kirkman et al. 2003, Freedman et al. 2001) are used when 
the CMB data are not included. We have also examined the case 
where the allowance for non-thermal pressure support has been 
doubled i.e. 1.0 < 7 < 1.2 (see text for details). 





Parameter 


Allowance 


Calibration/Modelling 


K 


1.0 ± 0.1 (Gaussian) 


Non-thermal pressure 


7 


1.0 < 7 < 1.1 


Gas depletion: norm. 


bo 


0.65 < b < 1.0 


Gas depletion: evol. 


a h 


-0.1 <a b < 0.1 


Stellar mass: norm. 


so 


0.16 ±0.048 (Gaussian) 


Stellar mass: evol. 


ct B 


-0.2 <a s < 0.2 


/gas(r ~ r 250 o) sl °P e 


>1 


0.214 ±0.022 (Gaussian) 


Standard prior Q b h 2 


Q h h 2 


0.0214 ± 0.0020 


Standard prior h 


h 


0.72 ± 0.08 


Weak prior Q,^h? 


n h h 2 


0.0214 ± 0.0060 


Weak prior h 


h 


0.72 ± 0.24 



Hubble parameter, 0.2 < h < 2.0. (Tests in which tensor 
components are included with ACDM models lead to similar 
results on dark energy, but take much longer to compute.) 



5 CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS 

5.1 Constraints on Q m from the low— z / gas data 

In the first case, we have used the Chandra / gaa data for 
only the six, lowest redshift clusters in the sample, with 
z < 0.15, to constrain the mean matter density of the Uni- 
verse. The restriction to low-z clusters minimizes correlated 
uncertainties associated with the nature of the dark energy 
component (dark energy has only a very small effect on 
the space-time metric over this redshift range; we employ 
a broad uniform prior such that 0.0 < Qa < 2.0) and ren- 
ders negligible uncertainties associated with the evolution 
of the depletion factor and stellar baryonic mass fraction 
(ah and a s ). Fig. [5] shows the marginalized constraints on 
Q m for a ACDM model with free curvature, using the stan- 
dard priors on Qhh 2 and h, for which we obtain a result of 
fi m = 0.28 ± 0.06. The full set of conservative systematic 
allowances, as described in Table [4] were included. 

The result on fi m from the six lowest redshift clusters is 
in good agreement with that obtained for the whole sample, 
as discussed below. It is also consistent with the result on 
Qm found from an analysis of all clusters except the six- 
lowest redshift systems, fi m = 0.29 ±0.06, i.e. the six-lowest 
redshift clusters do not dominate the fi m constraints. Note 
that the error bars on fl m are dominated by the widths of the 
priors on Qhh 2 and h and the magnitudes of the systematic 
allowances on K, b and 7, which are all at the ~ 10 — 20 
per cent level. In contrast, the statistical uncertainty in the 
normalization of the / gas (z) curve is small (Section 13. 1[) and 
has a negligible impact on the fi m results. 

The result on f2 m is consistent with previous findings 
based on / gas data (see references in Section 1) and indepen- 
dent constraints from the CMB (e.g. Spergel et al. 2007), 
galaxy redshift surveys (e.g. Eisenstein et al. 2005) and 




Figure 5. The marginalized constraints on O m from the Chandra 
/ gas data for the six lowest redshift clusters, using the non-flat 
ACDM model and standard priors on Q^h 2 and h. Uncertainties 
due to the evolution in b and s and the nature of the dark energy 
component are negligible in the analysis (although allowances for 
these uncertainties are included) . We obtain a marginalized result 
fim = 0.28 ± 0.06 (68 per cent confidence limits). 

other leading cosmological data. Note that the agreement 
in cosmological parameters determined from the / gas and 
CMB data argues against any unusual depletion of baryons 
within r25oo in hot, relaxed clusters (see e.g. the discussions 
in Ettori 2003, Afshordi et al. 2007 and McCarthy, Bower & 
Balogh 2007) 

5.2 Constraints on the ACDM model using the 
/ gas (+CMB+SNIa) data 

We next extended our analysis to measure f2 m and Qa for a 
non-flat ACDM model using the Chandra / gas data for the 
full sample of 42 clusters. The results are shown as the red 
contours in Fig. [5] Using the systematic allowances summa- 
rized in Table [4] and the standard priors on fit, ft 2 and h, 
we measure Q m = 0.27 ± 0.06 and Q A = 0.86 ± 0.19, (68 
per cent confidence limits) with % 2 = 41.5 for 40 degrees of 
freedom. The low x 2 value obtained is important and indi- 
cates that the model provides an acceptable description of 
the data (see Section f5. 31 below). The result on fl m is in ex- 
cellent agreement with that determined from the six lowest 
redshift clusters only ( Section [5TT} ■ The result is also consis- 
tent with the value reported by Allen et al. (2004) using the 
previous release of / gas data, although the more conserva- 
tive systematic allowances included here lead to the quoted 
uncertainties in Q m being larger by ~ 50 per cent. 

Fig. [7] shows the marginalized constraints on Qa ob- 
tained using both the standard and weak priors on Qhh 2 
and h. We see that using only the weak priors (Qhh 2 = 
0.0214 ± 0.0060, h = 0.72 ± 0.24), the / gas data provide a 
clear detection of the effects of dark energy on the expan- 
sion of the Universe, with Qa = 0.86 ± 0.21: a model with 
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Figure 6. The 68.3 and 95.4 per cent (1 and 2 a) confidence 
constraints in the fi m ,QA plane for the Chandra / gas data (red 
contours; standard priors on Q b h 2 and h are used). Also shown 
are the independent results obtained from CMB data (blue con- 
tours) using a weak, uniform prior on h (0.2 < h < 2), and SNIa 
data (green contours; the results for the Davis et al. 2007 compila- 
tion are shown). The inner, orange contours show the constraint 
obtained from all three data sets combined (no external priors 
on VL\>h? and h are used). A ACDM model is assumed, with the 
curvature included as a free parameter. 



Qa < is ruled out at ~ 99.98 per cent confidence. (Using 
the standard priors on Q^h 2 and h, a model with Qa < 
is ruled out at 99.99 per cent confidence; Table 0). The sig- 
nificance of the detection of dark energy in the / gas data 
is comparable to that of current SNIa studies (e.g. Riess et 
al. 2007; Wood-Vasey et al. 2007). The / gas data provide 
strong, independent evidence for cosmic acceleration. 

In contrast to the fi m constraints, the error budget for 
Qa includes significant contributions from both statistical 
and systematic sources. From the analysis of the full sample 
of 42 clusters using the standard priors on Q,bh 2 and h, we 
find J1a = 0.86 ± 0.19; the error bar is comprised of approx- 
imately ±0.15 statistical error and ±0.12 systematic uncer- 
tainty. Thus, whereas improved measurements of Q m from 
the / gas method will require additional information leading 
to tighter priors and systematic allowances, significant im- 
provements in the precision of the dark energy constraints 
should be possible simply by gathering more data (e.g. dou- 
bling the present / gaa data set). 

Fig. [6] also shows the constraints on f2 m and Ha ob- 
tained from the CMB (blue contours) and SNIa (green con- 
tours) data (Section I4.3[) . The agreement between the re- 
sults for the independent data sets is excellent and motivates 
a combined analysis. The inner, orange contours in Fig. [5] 
show the constraints on f2 m and Q.a obtained from the com- 
bined / gas ±CMB±SNIa data set. We obtain marginalized 
68 per cent confidence limits of fi m = 0.275 ± 0.033 and 



o <o 




Figure 7. The marginalized constraints on Qa determined from 
the Chandra / gas data using the non-flat ACDM model and stan- 
dard (solid curve) and weak (dashed curve) priors on fi^/i 2 and 
h. The / gas data provide a detection of the effects of dark energy 
at the ~ 99.99 per cent confidence level. 



fl a = 0.735 ± 0.023. Together, the / ga s±CMB±SNIa data 
also constrain the Universe to be close to geometrically flat: 
f2k = —0.010 ± 0.011. No external priors on Qbh and h are 
used in the analysis of the combined / ga s±CMB±SNIa data 
(see also Section [5. 6p . 

Finally, we have examined the effects of doubling the 
allowance for non-thermal pressure support in the clusters 
i.e. setting 1.0 < 7 < 1.2. For the analysis of the / gas data 
alone, this boosts the best-fit value of f2 m by ~ 5 per cent 
but leaves the results on dark energy unchanged. This can 
be understood by inspection of equation [3] and recalling that 
the constraint on il m is determined primarily from the nor- 
malization of the / g as curve, whereas the constraints on dark 
energy are driven by its shape (Section 14. 2[) . For the com- 
bined / g as±CMB±SNIa data set, doubling the width of the 
allowance on 7 has a negligible impact on the results, since 
in this case the value of f2 m is tightly constrained by the 
combination of data sets. 



5.3 Scatter in the / gas data 

Hydrodynamical simulations suggest that the intrinsic dis- 
persion in /gas measurements for the largest, dynamically 
relaxed galaxy clusters should be small. Nagai et al. (2007a) 
simulate and analyze mock X-ray observations of galaxy 
clusters (including cooling and feedback processes), employ- 
ing standard assumptions of spherical symmetry and hy- 
drostatic equilibrium and identifying relaxed systems based 
on X-ray morphology in a similar manner to that employed 
here. For relaxed clusters, these authors find that / gas mea- 
surements at r25oo are biased low by ~ 9 per cent, with 
the bias primarily due to non-thermal pressure support pro- 
vided by subsonic bulk motions in the intracluster gas. They 
measure an intrinsic dispersion in the / gas measurements of 
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Table 5. Summary of the constraints on cosmological parameters determined from the Chandra / gas data and complementary data sets. 
Error bars reflect the combined statistical and systematic uncertainties, incorporating the allowances and priors described in Section |4. 21 
For the low- z / gas data (z < 0.15), the constraint on f2 m is almost independent of the details of the dark energy component (Section [5TTJ. 
The SNIa(l) and SNIa(2) labels denote the supernovae samples of Davis et al. (2007) and Riess et al. (2007), respectively (Section 14.31 . 



Data 



Model 



Qb'i i h Priors 



COSMOLOGICAL CONSTRAINTS 

S~2de wo w c t 



low- 2 / gas 

/gas 

/gas 
/gas+CMB 

/ gas +CMB+SNIa(l) 



ACDM (0< Q A <2.0) standard 0.28 ± 0.06 

ACDM standard 0.27 ± 0.06 0.86 ±0.19 

ACDM weak 0.27 ± 0.09 0.86 ± 0.21 

ACDM none 0.28 ± 0.06 0.73 ± 0.04 

ACDM none 0.275 ± 0.033 0.735 ±0.023 



/gas 
/gas 
/gas + CMB 

/ gas +CMB+SNIa(l) 
/ gas ±CMB+SNIa(l) 



constant w (flat) 
constant w (flat) 
constant w (flat) 
constant w (flat) 
constant w 



standard 
weak 
none 
none 
none 



0.28 ±0.06 
0.29 ±0.09 
0.243 ±0.033 
0.253 ±0.021 
0.310 ± 0.052 



-1.14: 

-l.irt 



0.713 ±0.036 



,+0.27 
-0.35 
,+0.31 
-0.45 

-1.00 ±0.14 
-0.98 ±0.07 
-1.08±° ;» 



/ gas +CMB+SNIa(l) 
/ gas +CMB+SNIa(l) 
/ gas +CMB+SNIa(2) 



evolving w (flat) 

evolving w 
evolving w (flat) 



none 
none 
none 



0.254 ±0.022 

29+ 09 
0.287 ±0.026 



0.71 



+0.04 
0.05 



-1.05 
-1.15 
-1.19 



+0.31 
0.26 

+0.50 
-0.38 
-0.29 
-0.35 



-0.83.^ 

-,+0.70 
-1.30 
_ Q 33+0.18 
u - oo -0.34 



-0.80" 1 



~ 6 per cent, with an indication that the scatter may be 
even smaller for analyses limited to the hottest, relaxed sys- 
tems with kT > 5keV. Nagai et al. (2007a) also suggest 
that the true bias and scatter may be yet smaller if their 
simulations have underestimated the viscosity of the X-ray 
emitting gas01n contrast, for unrelaxed simulated clusters, 
Nagai et al. (2007a) find that / gas measurements are biased 
low by on average 27 per cent with an intrinsic dispersion of 
more than 50 per cent. Thus, the dispersion in / gas measure- 
ments for unrelaxed clusters is expected to be an order of 
magnitude larger than for relaxed systems. This is in agree- 
ment with the measurement of very low intrinsic systematic 
scatter in the / gas data for relaxed clusters reported here 
(see below) and the much larger scatter measured in previ- 
ous works that included no such restriction to relaxed clus- 
ters. Earlier, non-radiative simulations by Eke et al. (1998) 
also argued for a small intrinsic scatter in / gas , at the few 
per cent level, for large, relaxed clusters (see also Crain et al. 
2007). Likewise, Kay et al. (2004) measure a small intrinsic 
dispersion in / gas measurements from simulations including 
cooling and moderate star formation. 

The expectation of a small intrinsic dispersion in the 
/gas measurements for hot, dynamically relaxed clusters is 
strikingly confirmed by the present data. Even without in- 
cluding the allowances for systematic uncertainties associ- 
ated with 7, bo, Qb, s and a s described in Table [3] {i.e. 
keeping only the 10 per cent systematic uncertainty on the 
overall normalization, as described by K) the best-fitting 
non-flat ACDM model gives an acceptable x 2 — 41.9 for 40 
degrees of freedom, when fitting the full / gas sample. (The \ 2 
drops only to 41.5 with the full set of systematic allowances 
included; this small change in x 2 illustrates the degeneracies 

6 Recent work on the morphologies of X-ray cavities and Ha 
filaments suggest a relatively high gas viscosity (low Reynolds 
number) in nearby cluster cores (Fabian et al. 2003a, b, 2005, 
Ruszkowski, Briiggen & Begelman 2004, Reynolds et al. 2005). 



between the systematic allowances and model parameters.) 
The acceptable x 2 f° r the best-fitting model rules out the 
presence of significant intrinsic, systematic scatter in the 
current / gas data. This absence of systematic scatter is ob- 
served despite the fact that the root-mean-square scatter in 
the / gaa data is only 15 per cent. Moreover, the rms scatter 
is dominated by those measurements with large statistical 
uncertainties; the weighted mean scatter of the / gas data 
about the best-fit ACDM model is only 7.2 per cent, which 
corresponds to only 7.2/1.5 = 4.8 per cent in distance. 

5.4 Constraints on the constant w model using 
the / gas (+CMB+SNIa) data 

We have next examined the ability of our data to constrain 
the dark energy equation of state parameter, w. In the first 
case, we examined a geometrically flat model in which w is 
constant with time. Fig. [5] shows the constraints in the fl m , 
w plane for this model using the Chandra / ga s data and stan- 
dard priors/allowances (red contours), the CMB data (blue 
contours) and SNIa data (green contours). The different pa- 
rameter degeneracies in the data sets are clearly evident. 
For the / gaa data alone, we measure fi m = 0.28 ± 0.06 and 

The results for the three data sets shown in Fig. [S] 
are each, individually, consistent with the ACDM model 
(w — —1). The consistent nature of these constraints again 
motivates a combined analysis of the data, shown as the 
small, central (orange) contours. For the three data sets com- 
bined, we measure fi m = 0.253±0.021 and w = -0.98±0.07 
(68 per cent confidence limits). No priors on Q^h 2 and h are 
required or used in the combined / gas +CMB+SNIa analy- 
sis. The constraints on w from the combined data set are 
significantly tighter than 10 per cent. 

We note that our analysis accounts for the effects of 
dark energy perturbations, which must exist for dark en- 
ergy models other than ACDM; neglecting the effects of such 
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Figure 8. The 68.3 and 95.4 per cent (1 and 2a) confidence 
constraints in the Qm,w plane obtained from the analysis of the 
Chandra / gas data (red contours) using standard priors on Q^h 2 
and h. Also shown are the independent results obtained from 
CMB data (blue contours) using a weak, uniform prior on h (0.2 < 
h < 2.0) and SNIa data (green contours; Davis et al. 2007). The 
inner, orange contours show the constraint obtained from all three 
data sets combined: Q m = 0.253 ± 0.021 and w = -0.98 ± 0.07 
(68 per cent confidence limits). No external priors on Q^h 2 and h 
are used when the data sets are combined. A flat cosmology with 
a constant dark energy equation of state parameter ui is assumed. 

perturbations can lead to spuriously tight constraints (See 
Rapetti et al. 2005 for details). 

5.5 Constraints on the evolution of w from the 
combined / gas +CMB+SNIa data 

Fig. [9] shows the constraints on wo and w ct obtained from a 
combined analysis of / gas +CMB+SNIa data using the gen- 
eral, evolving dark energy model (Equation[6]l and assuming 
geometric flatness (fik = 0). The left and right panels show 
the results obtained for the two separate SNIa samples (Sec- 
tion 14. 3[) . Using the Davis et al. (2007) SNIa compilation 
(left panel), we find no evidence for evolution in the dark 
energy equation of state over the redshift range spanned by 
the data: the results on the dark energy equation of state at 
late and early times, wo = — 1.05^0 26 an d w ct = — 0.831°'^ 
(68 per cent confidence limits), are both consistent with 
a cosmological constant model (w — —1, constant). A 
similar conclusion is drawn by Davis et al. (2007) using 
SNIa+CMB+Baryon Acoustic Oscillation (BAO) data. 

We note, however, a hint of evolution in the dark energy 
equation of state when the Riess et al. (2007) 'gold' SNIa 
sample is used instead (right panel of Fig. . In this case, 
the marginalized constraints on dark energy at late and early 
times, as defined in Section l4~4l differ at the 2 — 3<r level. Sim- 
ilar indications are also apparent in the analysis of the same 



SNIa (+CMB+BAO) data by Riess et al. (2007). However, 
the analysis using the Davis et al. (2007) SNIa compilation 
(left panel), which includes the high-quality, high-redshift 
HST supernovae from Riess et al. (2007) and which shows 
no suggestion of a departure from the ACDM model, argues 
that the hint of evolution in the right panel of Fig. [9] may be 
systematic in origin (see also Riess et al. 2007 and Conley 
et al. 2007 for discussions). 

5.6 The degeneracy breaking power of the 
combined / gas +CMB(+SNIa) data 

The degeneracy breaking power of the combined / gas +CMB 
data set is evidenced in the left panel of Fig. 1101 which 
shows the constraints on fi m versus Ode for a ACDM model 
with free curvature for the CMB data alone (blue contours) 
and the combined / gas -|-CMB data set (orange contours). 
For the / gas +CMB data, we measure f2 m = 0.278±o;°5Q and 
= 0.732± 

0.046 (68 per cent confidence limits), with the 
curvature fik = — 0.011 jig '017 • As mentioned above, no ex- 
ternal priors on il^h 2 and h are required when the / gas and 
CMB data are combined. The degeneracy breaking power of 
other combinations of data with the CMB is discussed by 
Spergel et al. (2007). 

The right panel of Fig. HOI shows the constraints on the 
Hubble Constant, Ho, and mean baryon density, flbh 2 , de- 
termined using the flat, constant w model for the CMB data 
alone (blue contours) and the combined / gas +CMB data set 
(orange contours). The improvement in the constraints on 
these parameters determined from the / gas +CMB data over 
the CMB data alone is substantial. The tight constraints 
for the /gas+CMB data, H = 72.5 ± 4.6 kms" 1 Mpc -1 and 
fl^h 2 — 0.0223 ± 0.0007, demonstrate clearly why external 
priors on these two parameters are not required when the 
/ gaa and CMB data are combined. Indeed, the constraints on 
Ho and fl^h 2 obtained from the / gas +CMB data are signif- 
icantly tighter than the external priors on these parameters 
that are employed when the / gas data are used alone (Ta- 
ble |3I • Similar constraints on Ho and Fl^h 2 are presented 
by the WMAP team (Spergel et al. 2007) for flat ACDM 
models using various data combinations. 

Fig. I 111 shows the constraints on the dark energy equa- 
tion of state obtained from an analysis of the combined 
/gas+CMB+SNIa data set where the curvature is also in- 
cluded as a free parameter. The marginalized results for 
the constant w model (left panel), w = — 1. 08^0 19 an d 
fik = — 0-024^QQ^g, are comparable to those of Spergel et al. 
(2007; see their Fig. 17) from a combined analysis of CMB, 
SNIa and galaxy redshift survey data. The constraints for 
the non-flat evolving w model (right panel), though weaker 
than those for the flat model (Fig. [9}, remain interesting 
and are also consistent with a cosmological constant. As dis- 
cussed by Rapetti et al. (2005; see also Spergel et al. 2007), 
such results demonstrate the power of the / gaa +CMB+SNIa 
data to constrain the properties of dark energy without the 
need to assume that the Universe is flat. 

Using the non-flat evolving w model but fixing the tran- 
sition redshift Zt = 1 in Equation [6] we recover the model 
used by the Dark Energy Task Force (DETF) to assess the 
power of future dark energy experiments. The combination 
of current / gas +CMB+SNIa data provides a DETF figure 
of merit ~ 2. 
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Figure 9. The 68.3 and 95.4 per cent confidence limits in the (Q. rn \WQ,w e t) plane determined from the / gas +CMB+SNIa data using 
our most general dark energy model (Equation [5J with the transition scale factor marginalized over the range 0.5 < at < 0.95. The 
solid, purple contours show the results on (S7 m ,«>o)- The dashed, turquoise lines show the results on (Q m ,«! e t). The horizontal dotted 
line denotes the cosmological constant model (too = Wet = — 1). The left and right panels show the results obtained for the two SNIa 
samples: (Left panel) Davis et al. (2007) and (Right panel) Riess et al. (2007). A flat geometry (Q^ = 0) is assumed. The data provide 
no significant evidence for evolution in w and are consistent with the cosmological constant (ACDM) model [w = —1; Section |5.5|I . 
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Figure 10. The degeneracy-breaking power of the / gas +CMB data. Contours show the 68.3 and 95.4 per cent confidence limits determined 
from the CMB data alone (larger, blue contours) and combined / gas +CMB data (smaller, orange contours). (Left panel) the constraints 
on f2 m and Ode f° r the ACDM model with the curvature included as a free parameter. (Right panel) The tight constraints on Hq and 
£l^h 2 for the flat, constant w model, demonstrating why external priors on these two parameters are not required when the / ga a and 
CMB data arc combined. 



6 DISCUSSION 

The new Chandra / gas results and analysis presented here 
build upon those of Allen et al. (2004) and Rapetti et al. 
(2005). The present study includes 16 more objects, approx- 
imately twice as much Chandra data and extends the study 



beyond a redshift of 1. Our analysis includes a comprehen- 
sive and conservative treatment of systematic uncertainties 
(Section 14.21 see also Table [4}. Allowances for such uncer- 
tainties are easily incorporated into the MCMC analysis. 

As with SNIa studies, the / gas data constrain dark en- 
ergy via its effects on the distance-redshift relation to a 
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well-defined source population - in this case, the largest, dy- 
namically relaxed galaxy clusters - using measurements of 
a 'standard' astrophysical quantity - the ratio of baryonic- 
to-total mass in the clusters. Our results provide a clear 
and independent detection of the effects of dark energy on 
the expansion of the Universe at ~ 99.99% confidence for 
a standard non-flat ACDM model, an accuracy comparable 
to that obtained from current SNIa work (e.g. Astier et al. 
2006; Riess et al. 2007; Wood-Vasey et al. 2007; Miknaitis 
et al. 2007). Like SNIa studies, the / gas data trace the evo- 
lution of dark energy over the redshift range < Z < 1, 
where it grows to dominate the overall energy density of the 
Universe. Our results for the / gas data alone, and the com- 
bination of /gas+CMB+SNIa data, show that this growth is 
consistent with that expected for models in which the dark 
energy is a cosmological constant (w = —1). 

Despite some clear similarities, important complemen- 
tary differences between the / gaa and SNIa experiments ex- 
ist. In the first case, the physics of the astrophysical objects 
- large, relaxed galaxy clusters and SNIa - are very different; 
the fact that such similar cosmological results are obtained 
from the distance-redshift information for these separate 
source populations is reassuring. Future studies, combining 
the two techniques but using larger target samples, should 
open the possibility for precise distance-redshift measure- 
ments and good control of systematic uncertainties, employ- 
ing both kinematic and dynamical analyses (e.g. Rapetti et 
al. 2007; Riess et al. 2007 and references therein). 

An important strength of the / gas method is the tight 
constraint on fi m provided by the normalization of the 
/ gas curve; this breaks the degeneracy between the mean 
matter density and dark energy density inherent in the dis- 
tance measurements. Our result on Q m is consistent with a 
host of previous X-ray studies (Section 1). 

A further strength, which is of relevance when consider- 
ing observing strategies for future dark energy work, is the 
small intrinsic dispersion in the / gas distance measurements. 
SNIa studies have established the presence of a systematic 
scatter of ~ 7 per cent in distance measurements for indi- 
vidual SNIa using high quality data (Jha et al. 2007; see 
also e.g. Riess et al. 2004, 2007; Astier et al. 2006; Wood- 
Vasey et al. 2007). In contrast, systematic scatter remains 
undetected in the present Chandra / gas data for hot, relaxed 
clusters, despite the fact that the weighted mean statistical 
scatter in / gas data corresponds to only ~ 5 per cent in 
distance. This small systematic scatter for large, dynami- 
cally relaxed clusters (identified as relaxed on the basis of 
their X-ray morphologies) is consistent with the predictions 
from hydrodynamical simulations (e.g. Nagai et al. 2007a), 
although the results for both observed and simulated clus- 
ters are, at present, based on relatively small samples and 
more data are required. We stress that such small system- 
atic scatter is neither expected nor observed in studies where 
a restriction to morphologically relaxed clusters is not em- 
ployed e.g. compare the small scatter measured here with 
the much larger scatter observed in the studies of LaRoque 
et al. (2006) and Ettori et al. (2003); see also Nagai et al. 
(2007a). The restriction to the hottest, relaxed clusters, for 
which / gas is independent of temperature (Fig. [3J , also sim- 
plifies the determination of cosmological parameters. 

As mentioned above, the allowances for systematic un- 
certainties included in the analysis are relatively conserva- 



tive. Much progress is expected over the coming years in 
refining the ranges of these allowances, both observation- 
ally and through improved simulations. As discussed in Sec- 
tions 5.1 and 5.2, a reduction in the size of the required sys- 
tematic allowances will tighten the cosmological constraints. 
Improved numerical simulations of large samples of massive 
clusters, including a more complete treatment of star forma- 
tion and feedback physics that reproduces both the observed 
optical galaxy luminosity function and cluster X-ray prop- 
erties, will be of major importance. Progress in this area 
has been made (e.g. Bialek, Evrard & Mohr 2001, Muan- 
wong et al. 2002, Kay S. et al. , 2004, Kravtsov, Nagai & 
Vikhlinin 2005, Ettori et al. 2004, 2006, Rasia et al. 2006; 
Nagai et al. 2007a,b), though more work remains. In par- 
ticular, this work should improve the predictions for b(z). 
Further deep X-ray and optical observations of nearby clus- 
ters will provide better constraints on the viscosity of the 
cluster gas. Improved optical/near infrared observations of 
clusters should pin down the stellar mass fraction in galaxy 
clusters and its evolution. 

Ground and space-based gravitational lensing studies 
will provide important, independent constraints on the mass 
distributions in clusters; a large program using the Subaru 
telescope and Hubble Space Telescope is underway, as is sim- 
ilar work by other groups (e.g. Hoekstra 2007). Follow-up 
observations of the SZ effect will also provide additional, in- 
dependent constraining power in the measurement of cosmo- 
logical parameters (the combination of direct observations of 
the SZ effect using radio/sub- mm data and the prediction of 
this effect from X-ray data provides an additional constraint 
on absolute distances to the clusters e.g. Molnar et al. 2002, 
Schmidt, Allen & Fabian 2004; Bonamente et al. 2006 and 
references therein). Moreover, the independent constraints 
provided by the SZ observations should allow a reduction 
of the priors required in future work (e.g. Rapetti & Allen 
2007). 

In the near future, continuing programs of Chandra and 
XMM-Newton observations of known, X-ray luminous clus- 
ters should allow important progress to be made, both by 
expanding the / gas sample (e.g. Chandra snapshot obser- 
vations of the entire MACS sample; Ebeling et al. 2001, 
2007) and through deeper observations of the current target 
list. The advent of new, large area SZ surveys (e.g. Ruhl 
et al. 2004) will soon provide important new target lists 
of hot, X-ray luminous high redshift clusters. A new, large 
area X-ray survey such as that proposed by the Spectrum- 
RG/eROSITA project] could make a substantial contribu- 
tion, finding hundreds of suitable systems at high redshifts. 

Looking a decade ahead, the Constellation-X Observa- 
tory (Con-Xjf| and, later, XEUS0 offer the possibility to 
carry out precise studies of dark energy using the / gas tech- 
nique. As discussed by Rapetti & Allen (2007; see also 
Rapetti et al. 2006), the large collecting area and combined 
spatial/spectral resolving power of Con-X should permit 
precise / gas measurements with ~ 5 per cent accuracy for 
large samples (> 500) of hot, massive clusters (kT > 5keV) 
spanning the redshift range < z < 2 (typical redshift 



7 http:/ /www.mpe.mpg.de/projects.html#erosita 

8 http:/ /constellation. gsfc.nasa.gov/ 

9 http:/ /www.rssd.esa.int/indcx.php?project=XEUS 
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Figure 11. (Left panel) The 68.3 and 95.4 per cent confidence limits on the dark energy equation of state and curvature from the 
analysis of the / gas +CMB+SNIa data using the non-flat, constant w model. The SNIa compilation of Davis et al. (2007) has been used. 
The horizontal and vertical dotted lines denote the loci for cosmological constant models and geometric flatness, respectively, both of 
which are consistent with the data. (Right panel) The 68.3 and 95.4 per cent confidence limits in the (n m ;iAio,«J e t) plane determined 
from the / gas +CMB+SNIa data for the general dark energy model (Equation IB} with the curvature also included as a free parameter. 
Other details as in the left panel of Fig. [9] 



2 ~ 0.6). The predicted constraints on dark energy from 
such an experiment, assuming Planck priors (Albrecht et al. 
2006), have a DETF figure of merit > 20, which is compa- 
rable to other leading proposed dark energy techniques such 
as SNIa, cluster number counts, weak lensing and baryon 
acoustic oscillation studies. The high spectral resolution of- 
fered by the Con-X calorimeters will also permit precise 
measurements of bulk motions and viscosity in the cluster 
gas, addressing directly one of the main sources of system- 
atic uncertainty in the method. 

An ASCII table containing the redshift and / gas (z) data 
is available at http://xoc.stanford.edu or from the authors 
on request. The analysis code, in the form of a patch to 
CosmoMC, will be made available at a later date. 
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